DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Estimating  Surface  Currents  from  SAR 

Christopher  Wackerman 

1200  Joe  Hall  Drive,  P.O.  Box  990,  Ypsilanti  MI  48197 
phone:  (734)  480-5413  fax:  (734)  480-5367  email:  chris.wackerman@gd-ais.com 

Award  Number:  N00014-08-C-0561 


LONG-TERM  GOALS 

The  long-tenn  goal  of  this  program  is  to  be  able  to  utilize  SAR  imagery  to  estimate  surface  currents 
over  the  ocean.  With  the  availability  of  a  large  number  of  commercial  SAR  satellite  sensors,  it  may  be 
possible  to  have  an  operational  capability  to  generate  surface  current  maps  around  the  world's  oceans 
that  could  be  complimentary  to  the  operational  altimeters. 

OBJECTIVES 

There  are  two  objectives  in  this  program. 

(1)  Develop  a  capability  to  use  SAR  imagery  to  estimate  surface  currents. 

(2)  Develop  a  capability  to  assimilate  the  SAR  estimates  into  current  models  and/or  other 
observations  to  generate  a  wide-area  surface  current  estimate 

APPROACH 

(1)  Estimating  surface  currents  from  SAR 

Two  approaches  will  be  investigated  under  this  program.  One  will  be  the  inversion  of  radar  cross 
section  (RCS)  modulations  into  the  underlying  surface  current  gradient  that  has  caused  those 
modulations.  This  will  be  done  by  inverting  the  wave  action  balance  equation  that  determines  changes 
in  wave  spectra  due  to  interaction  with  surface  currents.  This  capability  has  been  generated  for  internal 
wave  signatures  by  another  ONR  program  (the  Non-Linear  Internal  Wave  Initiative).  Under  this 
program,  it  will  be  modified  for  larger-scale  current  patterns.  The  second  approach  will  be  to  generate 
a  two-dimensional  map  of  surface  wave  spectra  from  the  SAR,  and  use  the  spatial  changes  in  the 
spectra  to  estimate  surface  currents.  We  have  already  developed  tools  to  estimate  wave  spectra  under  a 
NOAA/NESDIS  program,  however  that  algorithm  will  be  modified  under  this  program  to  include  non¬ 
linear  transfer  functions.  We  anticipate  that  the  translation  of  wave  spectral  changes  into  surface 
currents  will  also  be  done  via  inversion  of  the  wave  action  balance  equation 

(2)  Assimilation  into  current  models 

It  is  not  clear  if  SAR  imagery  will  ever  be  able  to  provide  the  wide-area  view  of  surface  currents 
required  by  operational  users.  Thus  we  anticipate  that  we  may  need  to  assimilate  the  SAR-derived 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2Q  jq  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

4.  TITLE  AND  SUBTITLE 

Estimating  Surface  Currents  from  SAR 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROTECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

General  Dynamics  -AIS,1200  Joe  Hall  Drive,  P.O.  Box 

990, Ypsilanti, MI, 48197 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

4 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


currents  into  a  large  field-of-view  by  using  existing  models  of  current  fields.  We  propose  using  the 
DELFT3D  model  for  which  we  have  done  previous  inversion  work.  We  will  modify  this  previous 
work  so  that  we  can  detennine  the  model  inputs  required  to  reproduce  the  SAR-derived  currents,  then 
use  the  model  to  generate  the  large  field-of-view. 

WORK  COMPLETED 


We  have  finished  the  derivation  of  the  adjoint  expressions  for  the  wave  action  balance  equation  for 
estimating  the  two-dimensional  surface  current  field  from  observations  of  radar  cross  section  in  a  SAR 
image  and  we  have  implemented  it  in  software.  We  are  currently  testing  the  code  which  will  finish  the 
first  half  of  Objective  1 

RESULTS 


The  wave  action  balance  question  (shown  in  one  dimension  for  simplicity)  is 
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where  A  is  the  wave  action,  u  is  the  surface  current  in  the  x-direction,  cg  is  the  wave  group  velocity,  kx 
is  the  wave  number  in  the  x-direction,  and  Fs  represents  the  sources  and  sinks  of  the  wave  action 
(normally  a  wind  growth  tenn,  a  energy  dissipation  from  wave  breaking  term,  and  a  non-linear 
dissipation  term).  This  represents  the  changes  in  wave  spectra  across  a  scene  due  to  the  interaction  of 
the  waves  with  a  current  pattern.  We  have  developed  a  model  for  predicting  the  radar  cross  section 
(RCS)  of  a  SAR  image  from  the  map  of  2D  wave  spectra,  so  the  problem  to  solve  is  to  find  the  current 
pattern,  u(x),  such  that  the  resulting  wave  spectra  when  put  through  the  RCS  model  generate  the  “best 
“  fit  to  the  SAR  image. 

The  metric  we  will  minimize  is: 


J  =  JJJ(0(x)-M(x))2+D 
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where  D  is  a  Lagragian  multiply  to  make  sure  that  the  wave  spectra  satisfy  the  wave  action  balance 
equation  (this  is  become  the  adjoint  of  A),  M(x)  is  the  model  for  RCS  based  on  the  wave  spectrum  at 
location  x,  and  O(x)  are  the  observations  of  RCS  (i.e.  the  SAR  image).  We  want  to  find  u(x)  that 
minimizes  J,  and  we  want  to  do  that  via  a  search  operation.  This  means  that  we  need  to  know  the 
gradient  of  J  with  respect  to  u.  We  have  worked  through  the  variational  calculus  for  this  to  show  that: 
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where  D  must  satisfy  the  so-called  adjoint  equation 
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which  is  just  the  wave  action  balance  equation  with  a  different  right-hand  side.  Eq.  (4)  is  a  standard 
derivation  from  Eq.  (2);  Eq.  (3)  is  unique  to  the  surface  current  problem.  In  general  we  would  proceed 
by  starting  with  a  guess  for  u(x),  solving  Eq.  (1)  to  generate  the  wave  action,  A,  everywhere,  then 
solving  Eq.  (4)  to  generate  the  adjoint,  D,  everywhere,  then  plugging  A  and  D  into  Eq.  (3)  to  generate 
the  gradient  at  each  location  to  use  in  generating  the  next  guess  for  u(x).  However,  our  RCS  model, 
M(x),  is  a  two  scale  scattering  model,  which  means  that  there  is  a  Bragg  wave  scattering  component 
combined  with  a  modulation  of  the  Bragg  waves  by  the  long-scale  waves  that  is  dependent  on  integral 
properties  of  the  wave  spectrum.  This  means  that  the  term,  SM/SA  in  Eq.  (4)  is  dominated  by  values 
near  the  Bragg  wavelength.  This  is  illustrated  in  Figure  1  which  shows  an  image  of  3M/5A  (top) 
compared  to  an  image  of  A  (bottom).  This  means  to  first  order  we  only  need  to  solve  Eq.  (4)  for  these 
waves  (i.e.  for  the  spectral  locations  where  oM/SA  is  large).  These  waves  have  large  relaxation  rates 
and  thus  do  not  depend  significantly  on  the  transport  terms  in  the  wave  action  balance  equations.  We 
can  therefore  approximate  the  solution  to  Eq.  (4)  as 
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Note  that  this  is  the  exact  solution  for  a  Bragg  scattering  model  for  RCS  and  essentially  tells  us  how  to 
convert  errors  in  the  RCS  comparison,  0(x)-M(x),  into  gradient  estimates  (using  Eq.  (5)  put  into  Eq. 
(3)).  We  can  then  take  one  step  further,  and  note  that  the  long-wave  component  of  M(x)  is  dominated 
by  the  action  peak  (i.e.  the  bright  region  in  the  bottom  image  of  Figure  1)  whose  waves  have  small 
relaxation  rates  and  thus  just  spread  the  error  tenn  “upwave”  from  where  we  calculate  it  to  perturb  the 
estimate  of  surface  current. 


This  approach  has  been  implemented  and  is  currently  being  tested. 

IMPACT/APPLICATIONS 

If  successful,  the  resulting  algorithms  will  generate  an  operational  capability  to  estimate  surface 
currents  from  commercial  SAR  systems. 
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Figure  1:  Top  image  show  dM/dA  which  is  a  measure  for  how  much  a  change  in  one  spectral 
location  will  change  the  model  output.  This  term  is  dominated  by  the  region  around  the  Bragg  wave 
(bright  blobs  in  the  top  image)  since  the  model  is  dominated  by  Bragg  scattering.  Note  that  if  we 
only  used  Bragg  scattering,  the  top  image  would  be  two  dots.  Bottom  image  shows  the  action  (i.e. 

the  wave  spectrum)  which  is  dominated  by  the  wind  wave  peaks.  Thus  this  spectral  regions 
dominates  the  changes  in  the  model  due  to  changes  in  the  integral  measures  of  the  wave  spectrum. 


RELATED  PROJECTS 

The  proposed  research  will  start  with  algorithms  for  inverting  the  wave  action  balance  equation  to 
estimate  currents  that  have  been  developed  under  the  ONR  Non-Linear  Internal  Wave  Initiative.  The 
proposed  research  will  also  start  with  algorithms  for  assimilating  currents  into  DELFT3D  that  were 
developed  under  a  previous  ONR  program  entitled  "DELFT3D  Assimilation". 
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